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ABSTRACT 



We present evidence of magnetically- mediated disk accretion in Herbig Ae/Be 
stars. Magnetospheric accretion models of Balmer and sodium profiles calculated 
with appropriate stellar and rotational parameters are in qualitative agreement 
with the observed profiles of the Herbig Ae star UX Ori, and yield a mass ac- 
cretion rate of ~ 10~ 8 M yr _1 . If more recent indications of an extremely large 
rotation rate for this object are correct, the magnetic field geometry must devi- 
ate from that of a standard dipole in order to produce line emission consistent 
with observed flux levels. Models of the associated accretion shock qualitatively 
explain the observed distribution of excess fluxes in the Balmer discontinuity for 
a large ensemble of Herbig Ae/Be stars, and imply typically small mass accre- 
tion rates, < 1CU 7 M Q yr _1 . In order for accretion to proceed onto the star, 
significant amounts of gas must exist inside the dust destruction radius, which is 
potentially problematic for recently advocated scenarios of "puffed" inner dust 
wall geometries. However, our models of the inner gas disk show that for the typ- 
ical accretion rates we have derived, the gas should be generally optically thin, 
thus allowing direct stellar irradiation of the inner dust edge of the disk. 

Subject headings: accretion, accretion disks - - circumstellar matter — stars: 
emission-line — stars: pre-main sequence 



Introduction 



The class of Herbig Ae/Be (HAe/Be) stars was identified by Herbig (1960) in an attempt 
to find the analogues of T Tauri stars among objects of higher mass. This identification has 
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been supported by many studies over the last several years. Observations of mm-wave 
emission suggest the presence of circumstellar disks similar to those of the low-mass classical 
T Tauri stars (CTTSs) (Mannings, Koerner, & Sargent 1997; Mannings & Sargent 1997; Testi 
et al. 2000). Dusty disk models can in principle explain the infrared excess emission of many 
HAe/Be stars (Hillenbrand et al. 1992; Chiang et al. 2000; Natta et al. 2001; Dullemond, 
Dominik, & Natta 2001, DDN01). In addition, imaging observations in scattered light (e.g., 
Grady et al. 1999) also provide evidence for flattened or disk-like large scale dusty structures 
around a few HAe/Be stars. 

These studies naturally raise the question: are HAe/Be stars accreting from their disks? 
And if so, are the accretion rates high enough to modify stellar evolution, or to imply 
substantial migration of potential planet-forming material? The first attempt to derive 
HAe/Be accretion rates was made by Hillenbrand et al. (1992; HSVK92), who attempted to 
fit the observed infrared excess emission using steady accretion disk models and estimated 
much higher mass accretion rates (10~ 6 — 10 -5 M Q yr _1 ) than typical of CTTSs. However, 
HSVK92 also needed to invoke an inner disk hole, of order 0.1 AU in radius, to explain the 
decline in near-infrared excess emission shortward of A ~ 2 — 3/im. Hartmann, Kenyon, 
& Calvet (1993; HKC) argued that the accretion could not stop indefinitely at this inner 
radius, and showed that accretion at the rates inferred by HSVK92 would render any inner 
disk optically thick, eliminating the decline in infrared excess at short wavelengths. 

Natta et al. (2001) and DDN01 revisited this problem, this time developing models 
which explain the infrared excess emission of HAe/Be stars as the result of heating by 
irradiation from the central star, rather than local accretion energy release. These authors 
again assumed an inner disk hole to explain the decrease in near-infrared disk emission at 
short wavelengths, of size ~ 0.3 — 0.5 AU. The inner edge of the disk receives radiation 
from the star at near-normal incidence, and therefore becomes much hotter at a given radius 
than would be a geometrically thin, optically thick disk, irradiated obliquely. The increased 
heating thus causes the disk material at the inner edge to "puff up". DDN01 suggested 
that this expanded inner disk edge can quantitatively explain the magnitude of the near- 
infrared emission in HAe/Be spectral energy distributions (SEDs), without any accretion 
energy release. The location of the inner disk edge and its emission in this model is roughly 
consistent with recent interferometric observations (Millan-Gabet, Schloerb, & Traub 2001), 
which indicate near-infrared H and K band emission at radial distances well beyond that 
predicted by standard geometrically thin, optically thick, irradiated disks. 

While this model is attractive, the possibility of accreting material interior to the dust 
destruction radius must be addressed. All the low-mass T Tauri stars with near-infrared 
excesses exhibit accretion onto the central star, and it would be surprising if HAe/Be disks 
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with similar properties were completely inert. Indeed, heating of the inner disk edge by stellar 
radiation alone is likely to raise temperatures to the point that the fractional ionization is 
high enough that the magnetorotational instability (MRI) thought to drive disk accretion 
can operate (e.g., Balbus & Hawley 1998; Gammie 1996). In CTTSs, stellar magnetospheres 
truncate the inner disk, but as pointed out by HKC, the inner disk edge in the DDN01 model 
lies at such large radii that magnetospheric truncation is implausible, and in any event falls 
well outside of the corotation radius (where the disk Keplerian angular velocity matches 
the stellar angular velocity), which would therefore prevent accretion (e.g., Shu et al. 1994). 
Beyond this, there is direct evidence for ballistic infall in some HAe/Be stars (§2; Sorelli, 
Grinin, & Natta 1996; Natta, Grinin, & Tambovtseva 2000), which implies that material 
must extend well inward of the dust destruction radius. DDN01 argue that their model will 
hold even if accretion proceeds through the dust destruction radius, as long as the accreting 
gas is optically thin, but this assumption places severe constraints on accretion parameters 
(HKC) . Are accretion rates low enough that the inner accretion disk is optically thin? Or is 
there some other explanation of the near-infrared SED? 

In this paper we investigate accretion in HAe/Be systems and its implications for inner 
disk structure. We begin by analyzing a specific object of special interest, UX Ori, in terms 
of a magnetospheric accretion model similar to that developed by Muzerolle et al. (1998, 
2001). Next, we consider possible diagnostics of accretion from optical-ultraviolet continuum 
emission; these diagnostics limit accretion to rates much lower than those typically assumed 
by HSVK92. Finally, we investigate the inner disk structure in these systems. We show that 
for low accretion rates, the inner disk is likely to be optically thin, as assumed by the Natta 
et al. (2001) and DDN01 models; for somewhat higher accretion rates, the innermost disk 
may be optically thick but geometrically thin, allowing the direct irradiation of most of the 
"dust edge" of the disk. In this way we attempt to construct a more complete picture of 
HAe/Be systems. 



2. Accretion rates 

An important constraint on the optical depth in an accretion disk is the mass accretion 
rate M; higher M generally leads to greater optical depths in the disk (HKC). In this section 
we consider constraints on M from line profiles and UV/optical veiling, in the context of 
magnetospheric accretion. While early-type stars are thought to have much weaker magnetic 
fields due to dynamo activity than their cooler pre-main sequence counterparts, some main 
sequence A stars do exhibit ordered magnetic fields of strengths ~ kG (Mestel 1975; Borra, 
Landstreet, & Mestel 1982), which would be sufficient to truncate inner disks at a few stellar 
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radii (though not at the inner disk radii ~ 0.1 AU required by the DDN01 models; Hartmann 
1999). 

Moreover, magnetospheric accretion models are ideal for explaining high-velocity red- 
shifted absorption in some HAe/Be stars (Sorelli et al. 1996, and references therein). An 
alternate explanation of the ubiquitous redshifted absorption features seen in a number of 
HAe/Be stars is infall of cometary bodies (e.g., Grinin et al. 1994, G94; Sorelli et al. 1996; 
Grady et al. 1999). In this scenario, star-grazing planetesimals will evaporate, and, given a 
favorable viewing angle, redshifted absorption of neutral sodium and other atoms with low 
ionization potentials can result. This hypothesis provides a potential link to the older f3 
Pic objects with debris disks. However, a separate source, such as a wind or accretion flow 
is required to explain the line emission observed in stars such as UX Ori. Moreover, the 
infall signatures, while variable, are not episodic as in (5 Pic (e.g., Lagrange, Backman, & 
Artymowicz 2000, and references therein). Most significantly, Natta et al. (2000) found that 
the infalling gas is not strongly depleted in hydrogen, contrary to expectation if the source 
were evaporating comets (note that no circumstellar hydrogen absorption or emission has 
been observed in (3 Pic; Lagrange, Backman, & Artymowicz 2000). 

Below we consider magnetospheric accretion in more detail by modeling line profiles of 
one well-studied object, UX Ori. This star is particularly interesting for its large-amplitude 
optical variability, probably due to time-variable extinction by dust clouds (Grinin et al. 
1994; Natta & Whitney 2000). 

2.1. Line profiles 

UX Ori was observed with the echelle spectrograph on the 4m telescope at KPNO on 
September 17, 1999. The spectrum has a resolution of ~ 30, 000, and a wavelength coverage 
of 4300 - 9200 A, containing many features of interest such as the Balmer lines, He I, Na D, 
and the Ca II infrared triplet. The data were reduced using the standard IRAF echelle pack- 
ages. In order to provide pure emission profiles, the spectrum of standard HD 1280 (spectral 
type A2, v suu = 90 — 100 km s _1 ) was subtracted from the object spectrum, effectively re- 
moving photospheric absorption features (including the Balmer absorption wings). Figure 1 
shows the continuum-normalized profiles of nine atomic lines. The observed line profiles for 
UX Ori exhibit characteristics similar to those seen before (G94; Natta et al. 2000; Grinin 
et al. 2001). Few lines are in emission, but many exhibit redshifted absorption at supersonic 
velocities as high as several hundred km s -1 (e.g., Balmer, Na D, Ca II), indicating mass 
infall onto the star. 
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We constructed magnetospheric accretion models for comparison with the observed pro- 
files of UX Ori. The radiative transfer models are described in detail in Muzerolle et al. 
(2001). In summary, we assume a dipolar magnetic field geometry; gas accreting from a 
circumstellar disk falls ballistically along the field lines to the stellar surface. The density of 
the infalling gas is set by the geometry, velocity, and parameterized accretion rate. The inner 
and outer radius of the flow are free parameters (though the outer radius must be less than 
the disk corotation radius). Finally, we use a solid-body treatment for the rotation of the 
magnetosphere, including curvature of the field lines in the direction of rotation (Muzerolle 
et al. 2001). 

Model line profiles were calculated in the same manner as for the T Tauri comparisons 
with rotation, using the same gas temperature constraints found in Muzerolle et al. (2001) 
(i.e., in the range 6000-12,000 K, varying inversely with the density). Corresponding to the 
observed properties of UX Ori, an A2 star, the model stellar parameters were changed to 
M* = 3M Q and R* = 3R Q . We investigated two different values for the stellar rotation 
velocity, 70 and 140 km s" 1 , following v sini measurements from Bohm & Catala (1995) 
and Grinin et al. (2001), respectively, and given the expected edge-on orientation of the 
system. A more recent determination found vsini = 215 km s _1 (Mora et al. 2002); the 
reason for these hugely discrepant values is unclear. We discuss the effects of rotation and 
its implications on the size of the accretion flow below. 

Figure 2 shows a small sample of UX Ori model profiles, demonstrating the general 
effects of the accretion rate and inclination angle on the Balmer line profile shape and 
flux. Higher accretion rates produce broader emission wings from opacity broadening, as 
well as broader absorption from the increasing continuum opacity, which puts a limit on 
the peak flux such that the largest line fluxes occur at an intermediate accretion rate of ~ 
10~ 8 M yr _1 . Larger inclination angles result in more significant line asymmetries, including 
the appearance of broad, low-velocity redshifted absorption produced by the occultation of 
the hot star by the cooler outer parts of the magnetosphere. This particular feature, along 
with the stronger blue emission peak, is almost always seen in Ha: profiles of UX Ori and 
similar objects, and can only be explained by infalling material. 

The constraints our models can ultimately place on the accretion rate are limited by 
uncertain constraints on other model parameters, especially the gas temperature. For exam- 
ple, in some cases the same emission levels can be reproduced by increasing the temperature 
and decreasing the density (i.e., M), or vice- versa, beyond the T Tauri temperature con- 
straints. However, we find that M > 10~ 7 M yr _1 cannot reproduce the observations at any 
temperature: at T > 8000 K, the gas continuum optical depth becomes significant, result- 
ing in broad line absorption at Ha, while at lower temperatures, the Ha emission is much 
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weaker than observed. On the other hand, values of M much lower than 10~ 8 M yr" 1 also 
cannot reproduce the observed Ha emission, especially in the wings. In that case, simply 
increasing the gas temperature above that of the adopted constraints cannot result in much 
more line flux since the gas is already almost completely ionized. Ultimately, we believe a 
model-derived mass accretion rate estimate should be accurate to within about a factor of 
5. 

As already mentioned, another difficulty is the stellar rotation velocity, and the resulting 
limit on the size of the accretion flow. We find that the higher values of v sin % are problematic 
for the magnetospheric accretion model in that they result in very small corotation radii (2.1 
and 1.6 R*, respectively, for vsini = 140, 215 km s _1 ), and, hence, extremely small accretion 
flows since accretion cannot take place outside of corotation. For such small magnetospheric 
radii, there is not enough emission volume to match the observed Ha emission with any 
other combination of parameters; an example is shown in Figure 3. We obtain the best 
match to the observations using a larger magnetospheric radius, 3 i?*, roughly equivalent 
to the corotation radius for V* ~ 70 km s _1 . We note that the actual geometry of any 
magnetic fields that may be present in Herbig Ae stars is completely unknown. A "pinched" 
or nonaxisymmetric geometry that deviates from a standard dipole may allow for a model 
with a higher rotation velocity that has a larger emission volume and can reproduce the 
observations (Fig. 4 shows one possibility). Calculating the radiative transfer for such a 
complicated geometry, however, is beyond the scope of this paper. 

Figure 5 shows the observed and best model match profiles of Ha, H/3, and Na D for 
UX Ori. The magnetospheric infall models qualitatively match the observed line profiles in 
all three cases. The model accretion rate, 10 -8 Moyr^ 1 , is consistent with the upper limit 
we derive from the lack of hot continuum excess (see next section). The inclination angle for 
these models is 75°, in agreement with a nearly edge-on orientation of the system inferred 
from strong extinction events (e.g., G94). There is a central, slightly blueshifted absorption 
component in the Balmer lines that is not fit. This feature could be due to absorption from 
an accretion-driven wind exterior to the emission region (as seen in most CTTSs), which 
should be oriented roughly perpendicular to the accretion disk, resulting in small blueshifted 
absorption velocities given the edge-on orientation of the disk. The central absorption in the 
sodium lines is the interstellar component; there may also be some residual contamination 
from imperfectly subtracted night sky lines. The discrepancy in the level and velocity of 
the redshifted absorption components probably reflects deviations from the idealized flow 
geometry, as well as uncertainties in the gas temperature distribution. In any case, we feel 
the general agreement is reasonable, especially given the variable nature of the lines - note 
the similarity of these models to the multi-epoch spectra of UX Ori in G94. Similar Ha and 
Na D profiles have also been observed in spectra of BF Ori (de Winter et al. 1999), a Herbig 
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Ae star which exhibits photometric variability similar to that of UX Ori. Finally, we also 
note that our accretion rate estimate is comparable to that derived by Tambovtseva et al. 
(2001) for a model in which the infalling material lies essentially in the equatorial plane. 

As noted above, UX Ori is seen to to dim by up to several magnitudes in the optical 
(Herbst et al. 1983, 1994; Bibo & The 1990), accompanied by increased polarization (G94), 
which is interpreted as the result of an extinction event. During one such minimum, G94 
obtained an optical spectrum. The Ha profile exhibited a larger emission equivalent width 
(though with a lower total flux), with no absorption components, and a single emission peak 
blueshifted by nearly 50 kms" 1 . Meanwhile, the redshifted absorption in the Na D lines 
almost completely disappeared. These authors hypothesized that an optically thick dust 
cloud obscured most of the star during the minimum, while at the same time covering up 
much of the line emitting region, causing the change in profile shapes. Grinin & Tambovtseva 
(1995) calculated Ha profiles in rough agreement with the observations using a model of the 
circumstellar obscuration, and assuming gas in the inner disk with a velocity distribution 
similar to our magnetospheric accretion models. 

We attempted to simulate this obscuration with our models, to compare with both Ha 
and Na D. We included an opaque, nonemitting "screen" extending above and below the disk 
to a height of about 2.7 i?*, leaving just a small portion of the star and the magnetosphere 
unocculted. Such an occultation could be the result of density inhomogeneities or a warp 
in the disk or inner dust wall, enhanced by the nearly edge-on orientation of the system. 
The resulting profiles, shown in Figure 6, are in good agreement with the observations at 
photometric minimum. Ha has a much larger equivalent width, and shows a single peak 
blueshifted by about 50-100 kms -1 . The redshifted absorption in the sodium lines has 
completely disappeared. This behavior is easily explained as a result of the geometry of the 
accretion flow. The part of the flow moving away from the observer covers much of the star 
at the inclination angle used here (75°); the projection of this gas against the hot star results 
in redshifted absorption. The simulated obscuration occults most of this region, and the 
only part of the flow still observed is that moving towards the observer from behind the star. 
Thus, Ha shows a strong, blueshifted emission peak, and sodium does not show redshifted 
absorption (the emission measure of sodium is not high enough to produce any emission in 
this case). The similarity of the profiles in Figure 6 to those shown in Figures 4 and 5 of 
G94 lend further support for magnetospheric accretion taking place in UX Ori. 
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2.2. Magnetospheric accretion shock 

In the magnetospheric accretion model, material striking the stellar surface merges into 
the star through an accretion shock, from which most of the accretion luminosity is released. 
Models for the continuum emission from the accretion shock in CTTSs have been calculated 
by Calvet & Gullbring (1998, CG98). In these models, the gas in a plane-parallel accretion 
column, carrying energy flux JF, impacts the stellar surface and shocks to a temperature 
T s ~ 8.6 x 10 5 (M*/O.5M )/(R*/2R ) K, where M* and are the stellar mass and 
radius. The shock releases soft X-ray radiation which is absorbed by the accretion stream 
above and the stellar photosphere below the shock, producing optical and UV emission as it 
thermalizes. Approximately 3/4 of the emerging total column luminosity is emitted by the 
heated atmosphere below the shock, while the rest arises in the pre-shock region (CG98). 

CG98 and Gullbring et al. (2000) show that the excess continuum flux that veils the 
photospheric lines of T Tauri stars in the visible, and produces fluxes orders of magnitude 
above photospheric fluxes in the UV, can be explained in terms of the emission of accretion 
columns carrying energy fluxes of the order of T ~ 10 11 — 10 12 erg cm -2 s -1 , with a small 
surface coverage of accretion columns, / < 1%, where / is the filling factor. The high mass 
accretion rate stars, which show almost featureless spectra, could be explained with higher 
surface coverage, f < 10%. 

If HAe/Be stars undergo magnetospheric accretion, then an accretion shock must be 
formed at the stellar surface. Here, we calculate accretion shock models for stellar parameters 
appropriate to early A accreting stars, and with energy fluxes similar to those characterizing 
the lower mass accreting stars, following similar procedures to those in CG98. The total flux 
emerging from the star is 

F x = fFf + (1 - f)F? ot (1) 

where F phot is the flux from the undisturbed photosphere and F%° 1 = F^ p + F pres is the 
flux from the accretion columns, which in turn is the sum of the fluxes from the heated 
photosphere, F^ p , and the preshock region, F pres . The parameter / is the surface filling 
factor of accretion columns, estimated by the condition that the total luminosity carried 
by the accretion columns, L co i, is a fraction of the accretion luminosity. This fraction is 
determined by the disk truncation radius, Ri, through L co i = T x f'4nRl = (L acc , with 
L acc = GMM*/ R*, and ( = 1 — R*/Ri (CG98). In the calculations shown below, we take 
Ri = 2.5-R*, consistent with the magnetospheric infall models (section 2). 

Figure 7 shows the effect of the accretion shock on the spectrum of an underlying photo- 
sphere for an A2 star, with mass and radius M* = 3 M Q and R* = 3 R Q . Each row of panels 
corresponds to a different column energy flux, with T increasing from 10 10 ergcm~ 2 s _1 at 
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the top to 10 12 erg cm" 2 s~ x at the bottom. Each column corresponds to a different mass 
accretion rate onto the star, from 10~ 8 M yr _1 at the left, to M yr _1 at the right. 
For a given mass accretion rate and energy flux, the accretion columns must have surface 
filling factors indicated in each panel in Figure 7. Since L acc oc JF x /, for a given accretion 
luminosity, a lower energy flux results in a larger surface coverage. Since / cannot be > 1, 
high mass accretion rates cannot be achieved with low T columns alone. 

The contribution to the flux from the different zones in the accretion column, namely, 
the heated photosphere and the preshock region, are indicated in Figure 7, which also shows 
the emission from the undisturbed photosphere. The photospheric fluxes are taken from the 
stellar population library of Bruzual and Chariot (1993), while the fluxes from the heated 
photospheres are calculated here, assuming LTE and opacity sources as in CG98. This 
calculation does not include the hydrogen lines, which in a real spectrum pile up at the 
longward edge of the Balmer discontinuity, resulting in an effective shift of the wavelength 
of the jump, not present in the pure-continuum spectra of the heated photosphere. 

The main contribution to the flux in the visible and near-UV in Figure 7 comes from the 
optically thick heated photosphere below the shock. The radiative flux emerging per unit area 
from this region is F hp = F s + F* , where F s is the reflected irradiating flux from the shock 
regions above the optically thick photosphere, and F* the intrinsic stellar flux, entering the 
photosphere from below (CG98). In turn, F s ~ 3/4 T and F* = crT 4 ~ 3.6 x 10 11 (T./9000) 4 ; 
thus, if the energy fluxes carried by the accretion columns in HAe/Be stars are similar to 
those in CTTSs, T ~ 10 11 — 10 12 erg cm -2 s _1 , then the accretion energy input in the heated 
photosphere is comparable to or smaller than the intrinsic stellar flux. 

The similarity between the two fluxes entering the heated photosphere from above and 
below means that the effective temperature of the heated photosphere, T^ p = (-F^/cr) 1 / 4 , is 
only slightly higher than that of the undisturbed photosphere, as shown in Figure 7. Since / 
decreases as JF increases for a fixed mass accretion rate, the contribution to the emergent flux 
from the heated photosphere, fF hp (cf. eq.(l)), decreases as T increases. Note that even if 
the emission from the heated photosphere and the undisturbed photosphere are separated by 
some kind of "deveiling" procedure, the emission from the heated photosphere does not give 
a direct measurement of the accretion energy, except for very high values of T . This situation 
contrasts sharply from that in the cooler CTTSs, for which F* ~ 2 x 10 10 (T*/4000) 4 , and 
accretion dominates the column emission. 

However, there is a difference between the emission from the heated photosphere region 
and that from the undisturbed photosphere. The temperature gradient in regions with 
tross < 1 in the heated photosphere is shallower than that in the undisturbed photosphere, 
because of the extra heating from the shock above it. As a consequence, absorption features 
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formed in the heated photosphere will become weaker, and even disappear for high values 
of J 7 , in which the heated photosphere is nearly isothermal. So, for high M objects, with 
large /, the absorption lines will become increasingly veiled. The weakening of absorption 
features can be seen in Figure 7 comparing the Balmer jump (in absorption) in the heated 
photosphere for increasing values of T . The Balmer jump becomes weaker, and disappears 
for T ~ 10 12 ergcm -2 s _1 . As a result, the Balmer jump in the total emergent fluxes is 
increasingly filled in as T increases and the mass accretion rate, i.e. /, increases. 

This property suggests that the strength of the Balmer jump may provide a measure 
of the mass accretion rate in HAe/Be stars. Unfortunately, this is the least well-measured 
property of these stars. Garrison (1978, G78) carried out narrow (40 A) band photometry 
from 3400 to 8300 A for a sample of HAe/Be stars, and noted that the Balmer discontinuity 
was smaller than in main sequence stars of the same spectral type. G78 defined a measure- 
ment of the strength of the Balmer discontinuity, Db, as the difference in the magnitudes 
at both sides of the discontinuity, and a measurement of the excess emission, AD B , as the 
difference of Db for a given star to that of a standard of the same spectral type. Figure 8 
shows the distribution of values of ADb for stars with Db in Garrison's sample. We plot two 
histograms: one for the total sample of 16 stars 1 ; the other for the subsample of 12 stars 
with spectral types B5 - A5. The two distributions are essentially the same. Figure 8 also 
shows the predicted excess in the Balmer discontinuity for shock models with mass accre- 
tion rates between 10" 9 M yr _1 and 10~ 6 M yr _1 , and energy fluxes 10 10 erg cm -2 s~ 1 to 
10 12 erg cm -2 s -1 . We have calculated D B from the ratio of fluxes at 4000A and at 3640A, to 
make it consistent with Garrison's measurements. The Balmer jump excesses indicate low 
mass accretion rates < 10~ 7 M yr _1 for most of the stars in the sample. As stressed by 
G78, ADb is a reddening independent indicator, a property specially important for the "UX 
Ori" stars, affected by obscuration events in which it is difficult to disentangle absorbed and 
scattered light (cf. Natta & Whitney 2000). 

A few of the stars in Figure 8 have Balmer jump excesses consistent with high mass 
accretion rates, > 10 -6 M Q yr _1 . For these stars, G78 finds excess fluxes relative to the 
photosphere beyond A ~ 6500A. Our models predict excesses at long wavelengths, due to 
the optically thin emission of the pre-shock region, but the excesses are small. For instance, 
a model with M ~ 10~ 6 M Q yr _1 has ~ 0.1 mag excess in J — H and 0.15 mag excess in 
H — K. Such small excesses would be difficult to disentangle from observational errors and 
uncertainties due to spectral type mismatches, as well as excess emission from the disk. 

As discussed, spectra of HAe/Be stars with high mass accretion rates should have 



1 We have excluded Z CMa, which is an FU Orionis object (Hartmann & Kenyon 1996). 
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weaker, or "veiled", absorption lines than standard stars of the same spectral type. The 
abnormally weak features would be formed in the upper layers of the heated photosphere, 
and each line could be affected differently, depending on its exact height of formation. This is 
different from the case of low mass stars, where a "veiling continuum" affects all photospheric 
lines in a given wavelength interval, and can be easily extracted by subtraction of a standard 
spectrum from the target. For large M in HAe/Be stars, one would expect a large scattering 
in the values of the spectral types determined using different absorption features. Valenti, 
Johns-Krull, & Linsky (2000), in a study of the IUE short wavelength spectra of a sample 
of 74 HAe/Be stars, find that ~ 70% of the sample have spectral features in the ultraviolet 
consistent with those determined from optical spectra. In the case of UX Ori, discussed in 
section 2, the optical spectrum of UX Ori shows that only Ha is in emission, and there is no 
veiling in the absorption lines; analysis of the spectra as in Hernandez et al. (2002) yields a 
spectral type of A2 ± 1 subclass, using as indicators the Ca II K line, Fe II lines, Ca I lines, 
and H7. On the other hand, Grady et al. (1995) find that the 1560A photoionization edge in 
their IUE minimum light spectrum of UX Ori is consistent with spectral types A1-A2. The 
consistency of spectral type indicators from the optical to the ultraviolet stresses the lack of 
veiling in the absorption lines and indicates that the mass accretion rate is < 10~ 7 M yr _1 . 

In addition to continuum, the accretion column is expected to have and emission spec- 
trum of highly ionized metals (for instance, Si IV 1394, Si IV 1403, C IV 1549, C III 1909). 
These lines are observed in HAe/Be stars with spectral types later than A2, for which the 
weak photospheric emission below 1600 A does not obliterate the lines (Valenti et al. 2000); 
their strengths seem to be comparable to those of CTTSs. For hotter stars, these lines may 
appear in absorption, and even have a wind component observed at high resolution (Praderie 
et al. 1986; Bohm et al. 1996; Bouret et al. 1997; Grady et al. 1999). The wind may also be 
responsible for spectral features that appear deeper than expected from their spectral types 
in some regions of low resolution ultraviolet spectra (mostly corresponding to blends of Fe 
II lines, Valenti et al. 2000). In particular, we note that free- free emission from the wind 
may be responsible for the continuum excess between ~ 6500 A and ~ 1/iin in the stars 
with large Balmer jump excesses (Garrison 1978). Neither present shock models (see Fig. 7) 
nor wall emission (see below, section 3) can explain this excess. Appropiate modeling of the 
wind region is required to test this hypothesis further, and yield measurements of mass loss 
rates, thus providing lower limits to the mass accretion rates. 
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3. Accretion disk models 

Having determined that disks are accreting in HAe/Be stars, we now examine the struc- 
ture of the inner disks in order to address another question posed in the introduction, namely, 
whether or not they are optically thin, so that (1) they do not produce a large flux excess 
over the photosphere at near infrared wavelengths, and (2) they allow stellar radiation to 
reach the inner wall postulated by Natta et al (2001), Tuthill, Monnier & Danchi (2001) 
and DDN01 to explain the near-IR photometric excess and interferometric measurements. 
To illustrate the issues involved, we have calculated pure gaseous models for the inner disk 
regions of a typical HAe/Be system, with assumptions as those of DAlessio et al. (1999, 
2001). 

Optically-thin dust cannot be present interior to the dust destruction radius, Rd, due 
to the radiation field of the central object. In principle dust could exist if it were "shielded" 
from the stellar radiation by gas opacity, but in such a case the "irradiation surface" , defined 
as the surface at which the optical depth to the star (at characteristic frequencies) is ~ 1, 
has to be determined by the gas opacity. Thus any dust present at R < Rd would have to be 
below this surface. We have simplified the already complex treatment of this gaseous layer 
(see below) by assuming that dust is not present anywhere interior to the irradiation surface. 

If the sublimation temperature for silicates is taken as ~ 1500 K and the dust wall 
receives radiation frontally, as proposed by Natta et al. (2001) and DDN01, then the wall 
would be at a radius ~ 0.5 AU for the effective temperature and radius of a typical Ae 
star (see DDN01). The inner gaseous disk then corresponds to the region between the 
magnetospheric radius R ma g < 3 R* ~ 0.04 AU and ~ 0.5 AU. We calculate pure gaseous a 
irradiated accretion disk models for this region, using the typical HAe/Be stellar parameters. 

The vertical temperature structure is calculated using the formulation described by 
Calvet et al. (1991, 1992), modified to include viscous dissipation, which is quantified using 
the a prescription (Shakura & Sunyaev 1972). The stellar radiation field is treated as a 
parallel beam impinging on the disk surface at an angle cos" 1 /i relative to the normal. 
A fraction of the energy in the incident beam is scattered, maintaining the same frequency 
range as the stellar radiation field. The other fraction of the incident energy is absorbed 
by disk material, and we assume it is re-emitted at a frequency range characteristic of the 
local temperature. In the Calvet et al. (1991, 1992) formulation, the interaction between the 
direct and scattered stellar and disk radiation fields and the material is described by mean 
opacities, which we calculate self-consistently from monochromatic opacities using relevant 
sources for the conditions of the inner gaseous disk. 

For an HAe/Be star, approximately half of the stellar radiation is emitted at wavelengths 
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shorter than ~ 0.45 /zm, in the UV. In this range, we considered absorption due to H 2 (B- 
X, C + -X and C~-X, Abgrall et al. 1993) and CO (fourth positive system, Kurucz 1976) 
transitions, calculated with LTE molecular abundances. We found that the opacity due to 
the superposition of the rotational lines of CO and H 2 is so large for A < 0.2 \im that most 
of the stellar UV radiation is absorbed in the uppermost layers of the disk. A fraction of 
this UV radiation (with A < 0.11 /im) will photodissociate H 2 and CO (see for instance van 
Zadelhoff et al. 2003). Another fraction will be absorbed in CO and H 2 lines below the 
photodissociated layer. To calculate the structure of the gaseous inner disk, we make the 
crude approximation that the incident stellar radiation with A < 0.2 /im does not penetrate 
and heat the disk below a geometrically thin upper layer. In other words, we assume that 
all this radiation is absorbed and/or scattered by photodissociated molecules, or scattered 
in high-opacity lines that become saturated at the uppermost layers of the disk atmosphere. 

The remaining stellar radiation at A > 0.2 /im is mostly absorbed by TiO. We have 
incorporated the opacities of the a, (3, 7, 7', 5, e, and systems (Jorgensen 1994). We 
also included the opacity due to CN and CH (Kurucz 1995), although it turned out to 
be unimportant in the disk temperature range (~ 500 — 3000 K). At the characteristic 
disk wavelengths, the most important opacity sources are red and infrared bands of H 2 
(Auman 1967), pure-rotational bands of H 2 and OH (Tsuji 1966), and transitions of the 
electronic ground state of CO (Kirby-Docken & Liu 1978). We assume LTE and calculate 
the populations of TiO, H 2 0, C0 2 , C 2 , N 2 , 2 , CH, CN, CO, NH, NO, OH, H 2 , H^ , H~, H° 
using the method described by Mihalas (1967). The scattering of stellar radiation is mostly 
due to Rayleigh scattering of H 2 , with smaller contributions from Rayleigh scattering of H, 
He (Dalgarno 1962) and electronic scattering. 

Even with all these opacity sources, there is a lack of absorption in the wavelength range 
between ~ 0.2 — 0.45 /im. In this interval, radiation interacts with the gas only through 
Rayleigh scattering. Since scattering processes do not contribute to the disk heating, we 
make the additional approximation that only half of the incident stellar flux (mostly with 
A > 0.45 fiin or in the gaps between lines) contributes to the heating of the interior of the 
disk. 

The temperature as a function of the Rosseland mean optical depth r^, the mean scat- 
tering coefficient, and the ratio between the opacity to the stellar radiation and the Rosseland 
mean opacity, q, are calculated at each height and radius using equation (14) from Calvet et 
al. (1991), with a simultaneous solution of the hydrostatic equilibrium equation in the verti- 
cal direction. In order to evaluate the mean scattering coefficient and q (both assumed to be 
constant through the disk atmosphere), we computed the opacities at the upper temperature 
T and at a constant pressure P = 10~ 6 dyn/cm 2 , assuming both to be representative values 
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for the disk upper atmosphere. Given the proximity of the gaseous inner disk zone to the 
central star (R < 40 i?*), the curvature of the irradiation surface is not an important factor 
in the calculation of the stellar irradiation. So, we take the irradiation flux F irr and incident 
angle as those of a flat disk, 



the cosine of the angle between the incoming stellar beam and the normal to the disk. The 
factor 1/2 is introduced to account for our assumption that only half of the stellar radiation 
heats the regions of the disk below the uppermost layer (i.e., only half of the stellar radiation 
has A > 0.45/xm). 

With these assumptions and procedures, we calculate gaseous disk models for mass 
accretion rates, M = 1CT 9 , 10~ 8 , 1CT 7 and 1CT 6 M Q yr _1 , assuming a constant viscosity 
parameter a = 0.01. The resultant temperature structures of these models are shown in 
Figure 9. As shown by the dotted lines in Figure 9, the optically thin upper atmosphere has 
a temperature around ~ 2000 K, for all values of M. This is the equilibrium temperature 
corresponding to radiative heating by stellar radiation with A > 0.45 fim and radiative 
cooling, for a constant pressure of 10~ 6 dyn/cm 2 . This nearly constant temperature is due 
to the fact that in this temperature range for the assumed constant pressure, the cooling 
increases with temperature much faster than the heating, producing a thermostatic effect. 
Although the upper layer temperature is similar to or lower than the condensation limit for 
silicates, there cannot be dust in these (optically thin) regions because a dust grain would 
be heated by the stellar radiation to a temperature above its sublimation limit. As we have 
already mentioned, dust could exist at deeper layers if it were shielded by upper, optically 
thick gaseous layers; however, here we assume that the inner disk is free of dust at every 
height. 

The Rosseland mean optical depth of the gaseous inner disk increases with M. For 
M = 10 -9 M yr _1 , the disk is optically thin, and the midplane temperature is similar 
to or lower than the upper atmospheric temperature, as shown in Figure 9. At higher 
mass accretion rates, the innermost parts of the disk become increasingly optically thick. 
The radius inside which the disk is optically thick to its own radiation increases with M, 
reaching 0.25 AU for M = 10~ 6 MQyr^ 1 . As a result, the midplane temperature increases 
with M as shown in Figure 9. For M > 10~ 6 M Q yr _1 , it becomes higher than T ~ 10 4 K 
in the innermost regions, corresponding to the thermally unstable regime where hydrogen is 
partially ionized and the dominant opacity source is H~ (Bell & Lin 1994). 




(2) 



with 




(3) 
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Next, we calculate the height of the surface above the midplane, where the mean optical 
depth to stellar radiation is unity, z s . In our previous modeling efforts with dusty disks, 
we have used the Planck mean optical depth to estimate z s ; in the present case, the Planck 
mean optical depth to the stellar radiation (using the incidence angle for a perfectly flat 
disk) is larger than 10 3 , even for a low mass accretion rate like M = 10 _9 M Q yr _1 . However, 
the Planck mean opacity is not representative of the mean opacity for these cool gaseous 
disks; since the absorption of radiation is mostly due to a myriad of molecular rotational 
lines, stellar radiation is transported through the disk in the opacity gaps between lines. 
Thus, we have used a harmonic mean similar to the Rosseland mean, which gives a better 
representation of the actual transport of radiation. Figure 10 shows the height of the surface 
for models in Figure 9, calculated with the harmonic mean opacity to the stellar radiation. 
We find that for M = 10~ 9 M yr _1 , the harmonic mean optical depth to the stellar radiation 
is less than 1 for R > 0.1 AU. At R ~ 0.5 AU, zf s ~ 0.008, 0.02, 0.03 AU for M = 
10 -8 , 10~ 7 , 10~ 6 M yr _1 , respectively. This corresponds to z s /H(T c ) ~ 1, 2.6, 3.3, where 
H(T C ) is the scale height evaluated at the midplane temperature of each of these disk models 
(cf. Figure 9). 

At the dust condensation radius, the dust becomes the most important opacity source. 
For comparison with the gaseous disk, we estimate the height where the radial optical depth 
to the stellar radiation is one, behind the dusty wall. We assume that the gaseous inner 
disk does not absorb stellar radiation in the path between the star and the wall significantly 
(confirmed below). To calculate the optical depth to the stellar radiation in a disk annulus of 
radial thickness AR, we assume that the annulus is isothermal with a temperature T w ^1500 
K. Thus, the density can be approximated by a Gaussian p = p c exp(—z 2 /2H 2 ), where the 
midplane density is p c ~ S(i?) / \/2ttH(R), the surface density is S(i?) ps M /3nac s H(R), the 
scale height is H(R) = c s /Q, the sound speed is c s = a/ kT w / prriH , and Vt is the Keplerian 
angular velocity. The optical depth to the stellar radiation is 



'rad 



r+Ar 

p(r) x* dr, (4) 



where r is the radial distance from the star to a point (R, z) in the disk, such that dr = 
dR\/ R 2 + z 2 /R, and x* is the total mean dust opacity to the stellar radiation. Assuming 
that S, H, and are constant from R to R + Ai?, the integral can be approximated by 



r; ad ~x*<P>(R 2 + z 2 ) 1 ' 2 ^-. (5) 
The mean density < p > can be estimated as an average of the densities at (R, z) and at 
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(R + AR, z + Az), where Az w AR z/R, i.e., 



8ttH(R) 



e -z 2 /2H 2 _|_ e -z 2 (l+AR/R) 2 /2H 2 



(6) 



With AR/ R m 0.1, and estimating by assuming that the dusty wall contains grains with 
a maximum size a max = 0.1 /im, the heights where r* ad pa 1 are Zg USt /H(T w ) ~ 3.8, 4.3, 
4.8 and 5.3, corresponding to M — 10~ 9 , 10 -8 , 10~ 7 , 10~ 6 M yr _1 , respectively. These 
values are roughly consistent with the value found by DDN01, z^ ust ~ 5H(T W ), assuming an 
arbitrary power law surface density for the disk. 

The temperature of the dust wall T w is about 2 to 3 times higher than the temper- 
ature of the gas disk at the same radius (cf. Figure 9). Thus, z^ ust / z 9 s as ~ 10, 4, 3 
for M = 10~ 8 , 10~ 7 , 10~ 6 M yr _1 , as indicated in Figure 10. The fraction of the sur- 
face of the dust wall that receives direct stellar radiation is ~ 1, 0.9, 0.8, 0.7 for M = 
10~ 9 , 10~ 8 , 10~ 7 , 10~ 6 M yr _1 , and the portion of the wall in contact with the optically 
thick (to the stellar radiation) gaseous disk is small. Thus, at least for M < 10~ 7 M Q yr _1 , 
allowance for the stellar flux that cannot reach the wall because of absorption in the gaseous 
inner disk represents only a small correction to the Natta et al. (2001) and DDN01 treat- 
ments. 

The flux from the inner gaseous disk, in the most optimistic case of a pole-on orientation, 
is shown in Figure 11. We also include the contribution from the wall, which we have 
estimated with a blackbody of temperature 1500 K and an emitting area covering a solid 
angle 250 times the stellar solid angle (as expected for a vertical wall with height z s ~ 0.08 
AU, seen almost edge-on but without occultation of the emitting region). This blackbody 
contribution roughly matches typical HAe spectra (see HSVK92). The disk emission depends 
on M; in particular, the total emergent flux from the gaseous inner region in the if-band, 
relative to the star, is F K /F* tK ~ 0.5, 0.6, 1, 4 for M = 10~ 8 , 10~ 7 , 10~ 6 M yr" 1 . The inner 
gaseous zones of disks with M > 10~ 6 M yr" 1 are optically thick, and emit a large continuum 
excess in the near-IR, in agreement with Hartmann et al. (1993). For M < 10 -7 M yr _1 , 
the contribution to the total flux from the inner gaseous disk is negligible. The inner disk 
emits rotational- vibrational CO bands around 2.3 [im for M < 10 -7 M yr" 1 ; however, these 
CO emission bands are obliterated by the continuum emerging from the outer disk and the 
wall, in agreement with observed near-IR spectra of Herbig Ae/Be stars (unpublished IRTF 
SpeX data; Ishii et al. 2003). 
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4. Discussion 

We have shown that the magnetospheric accretion model can account for the overall 
line profiles of UX Ori and, by extension, similar systems. This requires substantial amounts 
of accreting gas which lie well within the radius at which dust is destroyed. Although we 
presume that the accreting gas departs from a disk structure within a few stellar radii of the 
photosphere, accounting for the large radial velocities observed, the co-rotation radius within 
which such magnetospheric accretion can occur is well interior to the dust destruction radius. 
Therefore, there must be a gaseous accretion disk between the dust opacity-dominated outer 
disk, responsible for most of the infrared excess emission, and the star (see Fig. 12). We have 
shown that for plausible parameters, this inner accretion disk can be optically thin. Even 
at somewhat higher accretion rates, an optically-thick inner disk can be geometrically thin, 
and thus not block most of the incident stellar radiation from reaching the outer dust disk 
This is consistent with the basic picture of Natta et al. (2000) and DDN01, in which there 
is a disk "wall" where dust sublimates. 

Our model, in which an accreting gas disk extends interior to the dust destruction 
radius, may help address or reduce stability problems. As pointed out by DDN01, the rapid 
increase in scale height of their model implies an inward radial gas pressure force at the dust 
destruction radius. DDN01 suggest that this might be countered by having the inner regions 
rotating faster than Keplerian, but this seems unlikely. The problem is alleviated to some 
extent in our model, in which there is an inner gas disk; the scale height variation is less 
than the photospheric surface z s variation. 

Our claim of magnetospheric accretion in UX Ori-type Herbig Ae/Be stars leads to the 
question of where the requisite magnetic activity originates. This is still far from clear, as 
the motivation for theoretical investigation has been lacking. Strong magnetic fields on the 
order of several kG have been measured in the lower-mass T Tauri stars (e.g. Johns-Krull, 
Valenti, & Koresko 1999). However, similar measurements for HAe/Be stars may prove 
exceedingly difficult, due to the significant rotational broadening of photospheric absorption 
lines. Given the parameters of our best line profile model, and using the relations given in 
Konigl (1991), we estimate a possible magnetic field strength of ~200 G (roughly half that 
if the small magnetospheric size is used) for UX Ori. We point out that at a younger age, 
further back along their radiative tracks, HAe stars once had cooler photospheres, and likely 
sustained convective dynamos (for example, intermediate-mass T Tauri stars such as T Tau 
and SU Aur, which have clear accretion signatures, may eventually become what we would 
classify as Herbig Ae stars). Thus, their current magnetic activity may be a relic of this 
earlier time. Detailed studies of the timescales of magnetic decay in such objects are needed 
in order to investigate this hypothesis further. 
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Fig. 1. — Line profiles from our echelle spectrum of UX Ori. Each profile was normalized by 
the continuum, and then subtracted by an A2 standard star spectrum. 
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Fig. 2. — Magnetospheric accretion model line profiles of Ha (solid line) and H/5 (dashed 
line) for the indicated inclination angles (in degrees, with zero being pole-on and 90 edge- 
on) and values of M (in units of M yr _1 ). All models were calculated with M = 3M Q , 
R = 3R , V* = 70 km s _1 , and magnetospheric radii R mag = 2.4 — 2.9. Maximum gas 
temperatures were T max = 8000 K for M = 10~ 7 M yr-\ 10,000 K for M = 10~ 8 M yr-\ 
and 12,000 K for M = 10" 9 Meyr" 1 . 
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Fig. 3. — Ha (solid line) and H/3 (dashed line) model line profiles for a smaller magnetosphere, 
with R mag = 1-35 — 1.85, corresponding to a stellar rotation rate of V* = 140 km s _1 . 
M = 10- 8 Moyr" 1 , T max = 10, OOOX. 



12 3 

R 

Fig. 4. — Comparison of different magnetosphere sizes (light and dark solid lines) used for 
the previous two figures, and an approximation of a "pinched" dipole configuration (dotted 
line), as discussed in the text. 
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Fig. 5. — Observed (solid lines) and model (dashed lines) line profiles for UX Ori. The 
model parameters are: M = f(T 8 M yr" 1 ; i = 75°; R mag = 2.4 - 2.9 R ; T rnax = 10, 000 K; 
Tphot = 8000 K; V rot = 70 kms" 1 . 
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Fig. 6. — Model profiles computed as in Figure 5 (dashed lines), compared to calculations 
with an opaque screen covering most of the star and emitting region (solid lines). 
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dM/dt=10" 




Fig. 7. — Emergent fluxes for A2 photospheres with mass accretion rates 10 -8 M Q yr _1 (left 
column), 10~ 7 M yr~ 1 (middle column), and 10~ 6 M yr _1 (right column). The energy flux 
of the accretion column is 10 10 erg cm -2 s -1 (upper row), 10 11 ergcm~ 2 s" 1 (middle row), and 
10 12 erg cm" 2 s -1 (bottom row). Filling factors for the accretion columns / are indicated in 
each panel. Heavy solid lines indicate total fluxes, light solid lines correspond to fluxes from 
the undisturbed photosphere, and dashed lines to the emission of the accretion column. The 
contributions from the shock regions are also indicated: heated photosphere (dot-dashed) 
and preshock (dotted). See §2.2. 
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Fig. 8. — Predicted excess in the Balmer discontinuity as a function of mass accretion 
rate, for different column energy fluxes: 10 10 ergcm^ 2 s _1 (dotted line), 10 11 ergcm _2 s _1 
(dashed line), and 10 12 erg cm -2 s" 1 (solid line). The distribution of observed excesses for a 
sample of 16 HAe/Be stars from Garrison (1978) is shown at the right. The shaded portion 
corresponds to the distribution of excesses for the HAe stars in the Garrison sample. Labels 
for the distribution are shown along the top axis. 
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Fig. 9. — Characteristic disk temperatures for mass accretion rates 10~ 9 , 10~ 8 , 10~ 7 , and 
1CT 6 M yr~ 1 , indicated in each panel: upper layer temperature T (dotted line), midplane 
temperature T mid (solid line), and viscous temperature T vis (dashed line). 
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Fig. 10. — Height of the inner gaseous disk surface for mass accretion rates 10 9 M yr 1 
(dash-dot), 10~ 8 M yr _1 (dot), 10~ 7 M^r" 1 (solid), 10~ 6 M^r" 1 (dash). The height of 
the dust wall is indicated by a solid square (see §3). 



31 



o 
o 



-9 



10 



-11 



-9 



10 



-11 



\\|r\ \ 
T \ - 

\ V"\ - 
\ \- 

: \ \: 

IdM/dt-10" 9 \ _ 

n. f i 1 i i i i i i i V 


: dM/dt-io- 8 \ : 


\ \ 

: \ : 

IdM/dt-10" 7 \ I 

~i i i i 1 i i i i i i i\ 


\ 

\ 

\ 

\ " 

IdM/dt-10" 6 \ : 

i i i i i i i i i i i\— 



0.5 1 0.5 1 
log A(/^,m) 



Fig. 11. — Emission from the inner gaseous disk (solid line) for mass accretion rates 
10 -9 , 10~ 6 , 10~ 7 , and 10~ 6 M Q yr _1 , indicated in each panel. The disks are seen pole-on 
and the fluxes are calculated at 440 pc. Emission by the dust wall, indicated by the dashed 
line, is shown as a blackbody at T = 1500K, covering a solid angle ~ 250 times the stellar 
solid angle, which is representative of observed spectra (see HSVK92). The SED of the 
stellar photosphere (dash-dotted line) was taken from the Bruzual & Chariot (1993) spectral 
library. The total flux is shown by the heavy solid line. The outer disk does not contribute 
at near-IR wavelengths (cf. Natta et al 2001). 
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Fig. 12. — A schematic view of disk accretion in HAeBe stars. Lines of constant gas density 
in the disk are indicated by thin curves; regions of relatively low (gaseous) optical depth 
are indicated by light shading, while optically-thick regions are shown by dark shading. At 
low accretion rates, the inner disk is optically thin (left panel), allowing direct irradiation 
of the dust "wall" (hatched shading). At higher accretion rates, even if the inner disk is 
optically thick (right panel) , it is geometrically thinner than the dust "wall" , again allowing 
for direct (normal incidence) irradiation. In this way the DDN01 model can be maintained 
while allowing for an inner disk which feeds magnetospheric accretion columns, funneling 
material down onto the central star (arrows). 



